Charm Coalescence at relativistic energies 
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The J/ip yield at midrapidity at the top RHIC (relativistic heavy ion collider) energy is calculated 
within the statistical coalescence model, which assumes charmonium formation at the late stage of 
the reaction from the charm quarks and antiquarks created earlier in hard parton collisions. The 
results are compared to the new PHENIX data and to predictions of the standard models, which 
assume formation of charmonia exclusively at the initial stage of the reaction and their subsequent 
suppression. Two versions of the suppression scenario are considered. One of them assumes gradual 
charmonium suppression by comovers, while the other one supposes that the suppression sets in 
abruptly due to quark-gluon plasma formation. Surprisingly, both versions give very similar results. 
In contrast, the statistical coalescence model predicts a few times larger J/ip yield in the most 
central collisions. 



A study of open and hidden charm production in 
nucleus-nucleus (A + A) collisions at RHIC (relativistic 
heavy ion collider) BNL is expected to shed light upon 
an important physical question of the space-time history 
of the charmonium formation. The standard "suppres- 
sion" approach is based on the idea of Matsui and Satz 
P|: charmonia are formed at the early stage of A + A 
reaction, the further evolution leads exclusively to their 
suppression due to interaction with initial nucleons from 
the colliding nuclei, secondary comoving hadrons, and/or 
deconfined medium. 

The idea of thermal J/ip production triggered the 
development of an alternative charmonium formation 
scenario, the statistical coalescence model (SCM) U0. 
Hidden charm mesons are assumed to be created at 
hadronization near the point of chemical freeze-out due 
to coalescence of charm quarks c and antiquarks c pro- 
duced at the initial stage. The distribution of c's and c's 
over different open and hidden charm species is given by 
the laws of equilibrium statistical mechanics. The SCM 
describes remarkably well the centrality dependence of 
the J/ip yield ^ as well as the transverse spectra |6J at 
SPS. 

A combination of the standard and SCM approaches 
as well as a nonthermal cc coalescence Q have been 
also considered. 

The preliminary RHIC data [9j on the J/ip rapidity 
density and its centrality dependence in Au+Au colli- 
sions at ^fs = 200 GeV from the PHENIX Collaboration 
have been already discussed in Refs. |lOUlll,ll2| . The final 
data ^3 , which became available recently, differ essen- 
tially from the preliminary ones. Because of low statis- 
tics, the data have rather large errorbars. The most likely 
value of J / ip yield in the most central collisions is not re- 
ported. Instead, only 90% confidence level upper limit 
is given. Still, the data appeared to be able to exclude 
the most extreme versions of the nonthermal coalescence 
scenarios 0. 

The aim of the present paper is to check whether SCM 



is tolerated by the data and compare its predictions to 
the standard suppression models. 

Let two nuclei A and B collide at impact parameter b. 
The number of produced J/ip mesons is given within the 
standard scenario by 0] 
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where CFj^ is the cross section of J/ip production in 
nucleon-nucleon (N + N) collisions, ?a(s) 1S the nuclear 
thickness function related to the nucleon density in the 
nucleus, and S(b, s) < 1 is a factor responsible for the 
J/ip suppression. 

At the very initial stage, charmonia experience absorp- 
tion, S = S abs , by sweeping nucleons of the colliding nu- 
clei (see, for instance, Refs. |l4 Il5|). Bound cc states 
are assumed to be absorbed in the so-called preresonance 
state, before the final hidden charm mesons are formed. 
The absorption cross section is therefore taken to be the 
same for all charmonia. The value a a b s = 4.4 mb [lfij 
follows from the most recent SPS data analysis and is 
close to the theoretical prediction of Ref. [rj|. We as- 
sume that the same value of a a bs prevails also at RHIC 
energies. 

Those charmonia that survive normal nuclear suppres- 
sion are subjected to the comover 0,0] or quark-gluon 
plasma (QGP) suppression [J^. Both suppression sce- 
narios describe successfully the centrality dependence of 
the J/ip yield in Pb+Pb collisions at the SPS. 

In the comover approach, an additional suppression 
factor appears: S = S abs S co 15j|. The suppression factor 
S co depends on the density of comovers and on an ef- 
fective cross section a co of J/ip dissociation by comovers 
(averaged over all comover species and all charmonium 
states contributing to the J/ip yield through their decays 
and also over particle momenta in the medium). The 
value a co = 0.65 mb [H from fits of new SPS NA50 data 
[2ll ] will be used in our analysis. We assume that the 
value of a rn remain the same also at RHIC. The charmo- 
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nium suppression at RHIC becomes, however, stronger, 
due to the higher comover density. 

There are two reasons for increasing the comover den- 
sity at RHIC relative to SPS. The multiplicity of pro- 
duced secondary hadrons per unit rapidity interval at 
midrapidity increases by a factor of about 1.5 from 
xfs = 17 GeV to y/s = 200 GeV already in elementary 
nucleon-nucleon collisions. Additionally, the deviations 
from the wounded nucleon model becomes stronger at 
higher energies. This increases the comover density in 
central nucleus- nucleus collisions. The centrality depen- 
dence of the number of light-flavored hadrons per unit 
pseudorapidity interval in Au+Au collisions at RHIC can 
be parametrized as [22j 
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where x = 0.11 for */a = 200 GeV |2J, N p (b) is the 
number of participants and N co u (b) is the number of col- 
lisions. Both are calculated in the Glauber approach. 

Calculating the centrality dependence of the J/ tp sup- 
pression, it is convenient to introduce an effective par- 
ticipant density in the plane transverse to the collision 



n*(b,s)= (1 — x)n p (b, s) + 2xn c (b, s) 



(3) 



Here n p (b,s) and n c (b,s) are, respectively, the densi- 
ties of nucleon participants and collisions in the trans- 
verse plane: N p (b) = J d 2 s n p (b,s) and N co u(b) — 
J d 2 s n co u(b,s). Note that the multiplicity of light- 
flavored hadrons J5J is proportional to N*(b) = 



comover density in the transverse plane, which is needed 
to calculate S co , is proportional to n*. 

The cross section of J/ip production per unit rapid- 
ity interval at midrapidity da^, 1 ^ / dy is the only free 

' y=o 

parameter of our fit. The PHENIX data on the J/ip 
multiplicity in p + p and Au+Au collisions pjjj are fit- 
ted simultaneously. The best fit, x 2 / n( V = 2.0 (3?|], is 
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is the branching ratio of J/ip decays into electron positron 
pair). The result is shown in Fig.^ 

The comover model has been historically referred to as 
a "hadronic" model. One might doubt whether the ex- 
trapolation of this model to the RHIC energies is legal. 
Indeed, the estimated energy density is extremely high 
even at SPS, so that hadrons can hardly preserve their 
individuality. The authors of the comover approach do 
not insist, however, on its hadronic interpretation (see, 
for instance We do not therefore make any as- 

sumptions about the nature of the comoving medium. 
We merely consider the comover model as an extreme 



scenario, which assumes a gradual increase of the char- 
monium suppression with growing energy density, with- 
out any abrupt changes of the absorptive properties of 
the medium. 

The QGP scenario of Ref.0] represents the opposite 
extreme: the charmonium suppression sets in, as soon as 
the energy density exceeds some threshold value. The 
excited charmonia, which contribute about 40% to the 
total J/ip yield, are suppressed at lower energy densities 
than directly produced J/ip's. We have updated the fit 
(l5| to the SPS data (new NA50 data jUf were added) 
using corrected values of the parameters, cr a fc s = 4.4 ±0.5 
mb and (jj T K / '(Tpy ~ 43.1, reported recently 0. Our 
results are 1%% = 2.99 fm~ 2 and n-i = 3.86 fm~ 2 . Here 
rxi (712) is the participant density in the transverse plane, 
corresponding to the threshold energy density at which 
excited charmonia (J/t/j's) are fully suppressed. 

Extrapolating to RHIC energies, one again has to take 
into account that the number of produced hadrons per 
unit rapidity and, consequently, the energy density of 
the produced medium grows with the collision energy 
and centrality. Due to the deviation from the wounded 
nucleon model J5J), the charmonium suppression sets in, 
when the effective participant density n*(b,s) © rather 

than the usual n p (b, s) exceeds the threshold values. The 
number of secondary hadrons per effective participant 
pair at y/s — 200 is higher than that at the SPS by a fac- 
tor of about 1.5. The critical energy density at RHIC is 
reached, therefore, at lower effective participant density: 
nl = ni/1.5 ps 2.0 fm" 2 and n* 2 = n 2 /1.5. w 2.6 fm~ 2 . 

Similarly as in the comover model, the J/ip production 
cross section is the only free parameter in the fit of the 
RHIC data. The minimum x 2 jridf = 2.2 is obtained at 



J d 2 s n* (6, s). Motivated by this fact, we assume that the B 
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= 5.0 x 10" 2 jib. The best fit of the 



QGP suppression scenario is also shown in Fig. 1. 

As was already noted above, the extrapolation of the 
standard suppression models from SPS to RHIC energies 
was based on the assumption that the value of the normal 
nuclear absorption cross section a a bs does not change es- 
sentially with the collision energy. Other viewpoints are 
also possible. The J / ip nuclear suppression mechanism at 
RHIC may be completely different from that at SPS [25| . 
This does not improve the agreement of the standard sup- 
pression models with the data, if the nuclear suppression 
becomes stronger at RHIC |26J. It is not excluded, how- 
ever, that the nuclear suppression may be even weaker 
|27|. A J/t/j measurements in d + Au collisions would 
clarify this point. 

Now we will check whether the statistical coalescence 
modelHH can be tolerated by the new data. In the SCM 
0, the total charm content of the final hadron system 
equals the number of c and c created at the initial stage 
of A + A reaction. Statistical laws control only the distri- 
bution of c and c among different hadron states in terms 
of the hadron gas (HG) model parameters: temperature 
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FIG. 1: The J/tp multiplicity per binary collision in the cen- 
tral unit rapidity interval at yfs = 200 GeV multiplied by 
the J/tp —> e + e~ branching ratio vs the number of nucleon 
participants. The most probable values (squares), statisti- 
cal errors (errorbars), and quadrature sum of statistical and 
systematic errors (boxes) are shown for the p + p and two 
centralities of Au+Au collisions. The 90 % confidence level 
(arrow) and its uncertainty (bracket) are shown for the most 
central collisions. 



T, baryonic chemical potential /i& and volume V. It ap- 
pears that the number of J ftp's produced by a statistical 
coalescence mechanism depends weakly on the thermo- 
dynamic hadronization parameters T and /xg. The J/tp 
yield is mainly defined by the average number of charmed 
quark-antiquark pairs A cg and by the hadronization vol- 
ume parameter V. If N cS is not much larger than unity, a 
proper account for the exact charm conservation becomes 
essential 0. This is crucial at SPS energies, where N cS is 
less than unity, and remains essential for moderate cen- 
trality in Au+Au collisions at RHIC. 

The SCM formula for the total (Air) J/tp yield that 
takes into account exact conservation of the number of cc 
pairs was obtained in Ref. 4] . In real experimental situa- 
tions, however, measurements are performed in a limited 
rapidity window. In the most simple case when the frac- 
tion of charmonia that fall into the relevant rapidity win- 
dow does not depend on the centrality, one can merely 
use the formula for the total yield multiplied by some fac- 
tor £ < 1. This approach was used [3 for studying the 
SPS data, where the multiplicity of light hadrons, which 
determine the freeze-out volume of the system, are ap- 
proximately proportional to the number of nucleon par- 



ticipants N p at all rapidities. At RHIC, the situation 
is different: the total (Air) multiplicity of light hadrons 
are approximately proportional to the number of partic- 
ipants, while at midrapidity, it grows faster [see Eq.J2J)]. 
The centrality dependence of charmonium production at 
different rapidities should, in this case, be also differ- 
ent. To compare the SCM prediction to the PHENIX 
data, which are related to the J/tp yield at midrapidity 
dNj/^p/dy, one has to derive a formula for the J/tp yield 
in a finite rapidity interval. 

Let £ak < 1 is the probability that a c quark, produced 
in a nucleus- nucleus collision, has rapidity within the in- 
terval Ay. The probability distribution of the number k c 
of c quarks inside the interval Ay for events with fixed 
total (47r) number N cS of cc pairs is given by the binomial 
law: 



f(k c \N c5 ) = 
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The probability distribution of the number k s of c's inside 
the interval Ay is assumed to be independent of k c |38| . 
It conforms to the same binomial law. Event-by-event 
fluctuations of the number of cc pairs N cS created at the 
early stage of A + A reaction in independent nucleon- 
nucleon collisions, are Poisson distributed: 



P{N C5 ; N cs ) = exp (~N c£ ) 



N r .,. ! 



(•5) 



The probability of cc coalescence is proportional to the 
product of their numbers and inversely proportional to 
the system volume. The proportionality coefficient de- 
pends on the thermal densities of the open and hidden 
charm, and is the same as in the case of the total char- 
monium yield Q. 

The average J/tp multiplicity at fixed values of k c and 
k e is therefore given by the formula |28| 
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[Deriving Eq.|0 we used the fact that the thermal num- 
ber of hadrons with hidden charm is much smaller than 
that with open charm.] Folding Eq.© with the binomial 
and Poisson distributions one gets 
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where no is the thermal density of all open charm 
hadrons and n j°L is the total thermal J / tp density (with 



decay contributions from the excited charmonium states 
included) . Both no and n j° 4 ^, are calculated in the grand 
canonical ensemble with the QGP hadronization param- 
eters T, ug, V A y found from fitting the data of light- 
flavored [3Jj hadron yields in the rapidity interval Ay. 
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The average number of cc pairs N cS is, however, related 
to their total (Att) yield. 

In Au+Au collisions at \fs = 200 GeV, the yield of 
light-flavored hadrons at midrapidity is fitted within the 
hadron gas model with T = 177 MeV and [i B = 29 MeV 
pjof . The centrality dependence of the volume is calcu- 
lated from 



Va v - 
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[the coefficient 1.2 is needed to recalculate the number of 
particles per unit pseudorapidity (77) interval to that per 
unit rapidity (y) interval [29j]. Here n c h is the charged 
hadron density calculated in the HG model. 

The average number of the initially produced cc pairs 
is proportional to the number of binary nucleon-nucleon 
collisions: N cS = N coU (b) crg N /ag^. The charm pro- 
duction cross section, <j^ 5 N , has been measured at RHIC 
by the PHENIX Collaboration j^- The result is consis- 
tent with PYTHIA calculations: aft N w 650 /xb. 

The SCM is applicable onlyto large systems: N p > 
100 in Pb+Pb at the SPS Uli Therefore, PHENIX's 
p+p point and the leftmost Au+Au point, corresponding 
to N p ~ 30, cannot be used in the SCM fit procedure. 
From this reason, we restrict ourselves only to a rough 
estimation of the SCM prediction for the J/ip yield at 
midrapidity at the top RHIC energy. 

We fix the charm production cross section in nucleon- 
nucleon collisions at its PYTHIA value, a^ s N = 650 /ib. 
There is no experimental data for the value of £ak=i, but 
one can roughly estimate it assuming approximately the 
same rapidity distribution of the open charm and J ftp's 
in p + p collisions. This leads to £,Ay=i ~ 0.3. The charm 
rapidity distribution in Au+Au collisions can be broader 
than in p + p due to rescattering of c and c by sweep- 
ing nucleons. This will not change our result essentially, 
however. The estimation of the total charm production 
cross section is based on the single electron measurement 
at midrapidity. Extrapolation to the total phase space 
has been done assuming that the charm rapidity distri- 
bution does not change from p+p to Au+Au. The charm 
production rate per binary collision at midrapidity was 
found to be independent of the centrality (at least within 
the present measurement accuracy) . This means that the 
value of the total charm production cross section would 
grow with the centrality, if there were a broadening of 
the rapidity distribution. Both effects, the decreasing 
of £Ay=i and the increasing of o~^ 5 N , nearly cancel each 
other in Eq. J7J and the prediction of SCM does not 
change significantly. 

The result is shown in Fig. ^ The SCM dependence 
of the J /if) rapidity density per binary collision on the 
centrality is almost flat at N p > 100, in contrast to the 
total J ftp yield, where a J/ip enhancement is expected 
[32]| . This difference appears because the hadronization 



volume at midrapidity, VA y =i, grows with the centrality 
faster than the total volume V. 

In conclusion, we have compared predictions of three 
different models for J/ ip production at the top RHIC en- 
ergy y/s = 200 GeV. None of the models are favored and 
none is excluded by the data. The statistical coalescence 
model predicts a few times larger J/ip production rate 
than the standard suppression models. Hopefully, mea- 
surements during the next Au+Au run will be able to 
clarify whether charmonia are formed only at the initial 
stage of the reaction (the standard suppression models) 
or production at the late stage via cc coalescence (SCM) 
is dominant. A crucial test for the SCM would be a mea- 
surement of the centrality dependence of tp' to J/tp ratio 
in Au+Au collisions. It should be constant (excluding 
the peripheral collision region) and equal to the value in 
equilibrium HG, if SCM is valid. 

The two standard charmonium suppression models, 
the gradual suppression by comovers and the abrupt sup- 
pression by QGP, give quite similar predictions. High 
quality data with small errorbars are needed to clarify, 
which of the models describes adequately the charmo- 
nium suppression process, if contribution from cc coales- 
cence is not significant. 

What is the charmonium production mechanism at 
SPS energies? Here both the standard suppression mod- 
els 0, ^3, 0| and SCM [j| , as well as their combination 
are also able to reproduce the data. However, SCM 
requires an essential (by the factor of about 3.5) enhance- 
ment of the open charm in Pb+Pb collisions. There is an 
indirect experimental evidence for such an enhancement 
[33]| and its possible mechanism has been considered [34| . 
Still, only a direct experimental verification [35| can give 
the final answer. 
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The J/ijj multiplicity 



in p + p is 



cross section by the standard formula 



related to the 



/dy 



y=0 

the total inelastic 



V<L do^/dy\ y=0 , where o& 
p + p cross section. 

Only statistical errors are taken into account in the cal- 
culation of x 2 - The most likely value of the J/ip yield 
in the most central collisions is not reported. Therefore, 
this point is not used in our fit procedure. 
This differs from Ref.|lO|. where exact equality, k c — ks, 
within the chosen interval Ay is assumed. In fact, net 
charm is exactly zero only in the total system. In any 
finite rapidity interval, event-by-event fluctuations with 
k c ^ k 5 are possible. 

At RHIC, the strangeness as well as all other conserved 
charges, excluding charm, can be safely considered in the 
grand canonical ensemble. 



